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Abstract 

Recently, observational searches for gravitational wave background (GWB) have been developed 
and given constraints on the energy density of GWB in a broad range of frequencies. These 
constraints have already resulted in the rejection of some theoretical models of relatively large GWB 
spectra. However, at 100 MHz, there is no strict upper limit from direct observation, though an 
indirect limit exists due to ^He abundance due to big-bang nucleosynthesis. In our previous paper, 
we investigated the detector designs that can effectively respond to GW at high frequencies, where 
the wavelength of GW is comparable to the size of a detector, and found that the configuration, 
a so-called synchronous-recycling interferometer is best at these sensitivity. In this paper, we 
investigated the optimal location of two synchronous-recycling interferometers and derived their 
cross-correlation sensitivity to GWB. We found that the sensitivity is nearly optimized and hardly 
changed if two coaligned detectors are located in a range ±0.2 m, and that the sensitivity achievable 
in an experiment is h^Q^Qg^ ~ 1.4 x 10^^. This would be far below compared with the constraint 
previously obtained in experiments. 
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I. INTRODUCTION 



There are many theoretical predictions of gravitational wave background (GWB) in a 
broad range of frequencies, 10^^^ — 10^° Hz. Some models in cosmology and particle physics 
predict relatively large stochastic GWB at ultra high frequency ~ 100 MHz; the quintessen- 
tial inflation model proposed by in which a large GWB spectrum is produ(^d_during 



the kinetic energy-dominated era after the inflationary expansion of the universe 



the violent reheating process a: 

n 

scenarios in string cosmology [9 



'ter inflation, a so-called preheating 

m 

101, 



MB, 



pre-big-bang 



black holes produced in the early universe 



ll[ |. the binary evolution and coalescence of primordial 
[isl . Q and their evaporation 15]. Recent 
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predictions of GW emission from black strings in the Randall-Sundrum model also gener- 
ate spectral features characteristic of the curvature of extra dimensions at high frequencies 



la, 113] . For the inquiry into high energy physics, testing these models with gravitational 



wave (GW) detectors for high frequencies is very important. 

Upper limits on GWB in wide-frequency ranges have been obtained from various obser- 



vations; cosmic microwave radiation at 10 



^18 



lO'i^Hz 



181 ]. pulsar timing at 10 ^ — 10 ^Hz 



|19IJ . Doppler tracking of the Cassini spacecraft at 10 ^ — 10 '^Hz 20], direct observation 
by LIGO at 10 - 10*H. H. ^He abundance due to big-baug uucleosyuthesis at frequeucies 

greater than 10 ^'^ Hz [22]. Nevertheless, as far as we know, no direct experiment has been 
done above 10^ Hz except for the experiment by A. M. Cruise and R. M. J. Ingley 23|]. They 
have used electromagnetic waveguides and obtained an upper limit on the amplitude of GW 
backgrounds, h < 10~^^ corresponding to h^Q^Qg^ < 10^^ at 100 MHz, where /iioo is the Hub- 
ble constant normalized with 100 kmsec"^ Mpc~^ and Qg„ is the energy density of GWB 
per logarithmic frequency bin normalized by the critical energy density of the universe, that 
is, 

^sw(/) = -p^. . (1) 
pcdlnf 

This constraint is much weaker than the constraints at other frequencies. Therefore, a much 
tighter boundary above 10^ Hz is needed to test various theoretical models. 

For our purpose of detecting GWB at 100 MHz, a detector with narrow bandwidth and 
high sensitivity is sufficient because we do not intend to know the exact waveform. In 
our previous paper [24], we considered three detector designs: a synchronous- recycling in- 



terferometer, a Fabry-Perot Michelson interferometer and an L-shaped cavity Michelson 
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interferometer, and investigated their GW responses. As a result, we found that the syn- 



chronous recychng interferometer (SRI) 25| is the most sensitive at lOOMHz, in which the 
angular mean sensitivity is slightly worse than that in low frequencies because the detector 
size is comparable to GW wavelength. 

To detect GWB with smaller amplitude than detector noise, one has to correlate sig- 
nals from two detectors in order to distinguish the GW signal. The analytical method has 



been well developed by several authors 



26 



21, 



28j. In these references, they assume that 



GW wavelength is much larger than the detector size, which is a so-called long-wave ap- 
proximation. However, it is not valid in our situation around 100 MHz, in which the GW 
wavelength is comparable to the detector size. This means the relative location of the two 
SRIs significantly affects the correlation sensitivity to GWB and the response function of 
a detector should be taken into account properly. Thus, in this paper, we will extend the 
previous analytical method of correlation, including the response functions of the detectors, 
and investigate the dependence of the sensitivity on the relative location of two detectors. 

This paper is organized as follows. In Sec. II, we will briefly review the correlation 
analysis for GWB and extend it including a full detector-response function. In Sec. Ill, we 
will explain a synchronous recycling interferometer and its response function. Sec. IV is the 
main section of this paper and is devoted to the evaluation of an overlap reduction function 
and the consideration of dependence of sensitivity on the locations of detectors. We will 
calculate the sensitivity of two SRIs to GWB in Sec. V and conclude in Sec. VI. 

II. CORRELATION ANALYSIS 

In this section, we will briefly review the formalism of correlation analysis for GWB 



221 . |26| | and redefine a part of the formalism, including detector response functions, which 
are necessary to take the finite size of detectors into account. 

Let us consider the outputs of a detector, s{t) = h{t) + n{t), where h{t) and n{t) are 
the GW signal and the noise of a detector. At generic point X, the gravitational metric 
perturbations in the transverse traceless gauge are given by 

h(t, ±) = J2 dn df h,{f, O) e^-^^*-"-*/'^) e,(0) , (2) 

pJs^ J~oo 

where c is speed of light, O is a unit vector directed at GW propagation and hp{f, ft) is 



the Fourier transform of GW amphtude with polarizations p = +,x. Polarizarion tensors 
ep(ri), can be written as 

e+(fi) = m(8)m — n(8)n, 
ex(r2) = m(g)n + n(g)m, 

where the unit vectors rh, fi are orthogonal to O and to each other. 

In this paper, we assume that GWB is (i) isotropic, (ii) unpolarized, (iii) stationary, and 
(iv) Gaussian, and (v) has small amplitude compared with that of noise, \h{t)\ <^ 1^(^)1- 



22|. 



These assumptions (i) - (iv) are discussed in [261], and (v) is expected at high frequencies 
around 100 MHz because there exists an indirect upper limit on GWB by big-bang nucle- 
osynthesis. These assumptions (i) - (iv) are expressed by 

{h;{f, n)h,if', n')) ^ 6{f - f')^s\n, il') V^'5,(/), (3) 

where 6'^{Q,Q') = 5{(j) — (f)')6{cos6 — cos6'), Sh{f) is the one-sided power spectral density 
and is defined by Eq. ([3]), and (■ ■ ■ ) denotes ensemble average. The power spectral density 
Sh{f) is related to Qg^j by 

^.M) = (^] fs.if) , (4) 

where Hq is the Hubble constant 

GW signal h{t) from a detector is given by D(/, fl) : h(/, fl), where the symbol : denotes 
contraction between tensors, and D(/, fi) is a so-called detector tensor, which describes the 
total response of a detector and maps the gravitational metric perturbation to the GW 
signal from a detector. We define it including detector response functions as 

D(/, n) = ^ [(u ® u)T(/, rj ■ u) - (v ® v)T(/, O • v)] . (5) 

Here u and v are unit vectors. We assume that they are orthogonal to each other and are 
directed to each detector arm. The function T is a detector response function, which we 
will introduce below, that describes the effect of finite arm length on propagating light and 
gives unity in low frequency limit. In the detector whose arm length is much smaller than 
the wavelength of GW, this function is approximated to unity, while in our detector whose 
detector size is comparable to GW wavelength, the function significantly affects the response 
of the detector. Using Eqs. ([2]) and (j5]), GW signal h{t) can be written as 

hit, x) = y2 dn df h,if, n) e2-/(*-"*/^)F,(/, n) , (6) 



where an angular pattern function of a detector Fp{f, f2) is defined by 



Fp{f,fl) = -D{f,n):e,{a). 



(7) 



Cross-correlation signal S between two detectors is defined as 




■T/2 



■T/2 



dt' si{t)s2it')Q{t-t'), 



(8) 



where Si and S2 are an output from each detector, T is observation time. Q{t — t') is an 
arbitrary real function, which is called an optimal filter. We will determine its form below 
so that signal-to- noise ratio (SNR) is maximized. Fourier transforming si{t) and S2{t), one 
can obtain 



tively. 5t(/) is the finite-time approximation to the Dirac delta function defined by 



In the above derivation, we took the limit of large T for one of the integrals. This is justified 
by the fact that, in general, Q{t — t') rapidly decreases for large \t — t'\. The correlation 
signal obtained above ideally has a contribution from only the GW signal since we assume 
that noise has no correlation between two detectors. Thus, we take ensemble average of Eq. 
dn]) and obtain the signal from GWB, 





where S2{f) and Q{f) are the Fourier transforms of Si{t), S2{t) and Q{t — t'), respec- 





(10) 



Substituting the Fourier transform of Eq. Qj, 




(11) 



into Eq. ( ITOi) . and using Eqs. ([3]) and (jll), one can obtain 




(12) 



Here we defined the overlap reduction function. 



7 



p 



(13) 



6 



— * — * — * _ 

where the separation of two detectors is AX = Xi — X2. The factor in front of Eq. fll3p is 
a normahzation and the overlap reduction function gives unity in low frequency limit. This 
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In the references, the overlap 



definition is slightly different from that in other papers 
reduction function is defined as meaning how GW signals in two detectors are correlated, 
and equals unity for colocated and coaligned detectors, while our definition includes the 
detector response functions T, as one can see explicitly in Eq. ([5]), and does not give unity 
even for colocated and coaligned detector at high frequencies. Namely, we regard the loss 
of GW signals of detectors as the reduction of overlap between two detectors. We will see 
below that this is important for the detector whose size is comparable to GW wavelength. 

Next, we will calculate the variance of a correlation signal. In this paper, we assume that 
noises in two detectors do not correlate at all and that the magnitude of GW signal is much 
smaller than that of noise. Consequently, the variance of correlation signal is 

Then, using Eq. ([9l), it follows 

/OO POO 
df / dfQUWW) {W)hU')) {hUKU')) 
-co J —00 

rji poo 

^7/ dfp.iifmmiQUii', (14) 
^ J —00 

where the one-sided power spectrum density of noise is defined by 

Now let us determine the form of the optimal filter Q{f). Equations f|T2l) and f|T^ are 
expressed more simply, using an inner product 

/OO 
dfA*if)Bimi\f\)p,i\f\), 
-00 

as 

^ 207r2n^'l/PA(|/|)P2(|/|)y' ' ^ ^ 



4 

From Eqs. f|T5l) and f|T6l) . SNR for GWB is defined as SNR = yu/o". Therefore, the optimal 
filter function turns out to be 



Qif) = A 



7(/)^gw(|/|) 

ifi'PiiifmnfW 



with an arbitrary normalization factor A. Applying this optimal filter to the above equations, 
we obtain maximal SNR 



SNR 



3m 



.1/2 



107r2 

3.19 X 10"^^ 



df 



In the transformation from the first to the second, we used Hq = 100 /iioo km s~^ Mpc ^ 



(17) 



3.24 /linn X 10 



-18 ^-1 



III. SYNCHRONOUS-RECYCLING INTERFEROMETER 

In this section, we will describe the detector response function T(/, f2) and spectrum 
density of noise Pi, i = 1,2 of a synchronous-recycling interferometer (SRI), which are 
needed to calculate the sensitivity to GWB. The design of SRI is shown in Fig. [H which was 



first proposed by R. W. P. Drever in 25|] . Laser light is split at a beam splitter and is sent 
into a synchronous-recycling cavity through a recycling mirror, which is the mirror Mi in Fig. 
[H The beams circulating clockwise and counterclockwise in the cavity experience GWs and 
mirror displacements, leave the cavity, and are recombined at the beam splitter. Then, the 
differential signal is detected at a photodetector. The advantage of SRI is that GW signals 
at certain frequencies are accumulated and amplified because the light beams experience 
GWs with the same sign of phases during round trips in the folded cavity. Consider GW 
propagating normally to the detector plane with an optimal polarization. In this case, the 
GW signal is amplified at the frequencies / = {2n — 1) x c/AL, n = 1, 2, ■ ■ ■ , where L is 
the arm length [30!]. The disadvantage of SRI is less sensitivity for GWs at low frequencies, 
/ < c/4L, because the GW signal is integrated in the cavity and canceled out as the 
frequency decreases. 



In our previous paper [2j], we derived the GW response of an SRI and found that it 
can be written in the form 5$(/, Q,) = a{f) S(j){f, O), where 50(/, tl) denotes the Fourier 
component of phase shift due to GW during the round trip of light in a recycling cavity and 
a{f) denotes an optical amplification factor of light in the cavity. The concrete expressions 
are written as 

"^^^ " ~ {Rf - RE)il - RfRe e-^-'f^) ' 
8 
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FIG. 1: Design of synchronous recycling interferometer (SRI). 
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FIG. 2: (color online) Arm response functions as a function of /r. Each plot is for arccos(r2 • u) 
7r/2, 7r/3, 7r/6, respectively, as shown by plot labels in the figure. 



'l-e 



-4mfT\ 



2nf 



V (g) V 



1- (n- v)2 

u ® u 



{ sin(27r/r) -i(i7 ■ v)(e 



-2lTifT fl-v 



cos(27r/r))} 



{ sin(27r/r) -i(i7-u)(e 



—2nifT tin 



-cos(27r/r))} 



, (19) 



where t = L/c, uj is the angular frequency of laser, X is a position vector of the mirror 
Ml, Rp is the amplitude reflectivity of a front mirror, and Re is the composite amplitude 
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reflectivity of tlie otlier tliree mirrors in tlie cavity. Tp is given by T|. = 1 — since we 
are assuming none of tlie mirrors have losses. The unit vectors u, v and tl are the same as 
those in Sec. [TTl Here we define an arm response function 



r(/, n ■ u) 



27r/r{l-(0-u)2} 



sin(27r/r) -i{n- u){e-2^^^""-" - cos(27r/r)} 



(20) 



so that T gives unity in low frequency limit, as shown in Fig. [2j Using this response function, 
we can rewrite Eq. (I19p in a simple form. 



50(/,O) = 2 cure 



-27ri/A-X/c 



1 -e" 



-47r^/T^ 



(u ® u)T(/, n■u)-{^r(E) v)r(/, U ■ v) 



Therefore, comparing Eq. (!2T|) with Eq. (fTTl) . we obtain 

= K{f)a{f)h{f). 



(21) 



(22) 



where k(/) = 2co'r(l — e ^'^^Z^). To identify /i(/) and 5$ (/), we incorporate the extra factor 
ft(/) a{f ) into the noise spectrum|3l[|. that is. 



/?W7) = l'^(/)«'(/)rSMnoise(/) 



'2ncu 



(23) 



where h is the reduced Planck costant, rj is the quantum efficiency of photodetector, Iq is 
laser power, and Pqnoise is the noise due to vacuum fluctuations. We assumed that shot noise 
is a dominant noise source around 100 MHz. In Eq. fl23l) . «(/) is replaced with 



«'(/) 



ReTp 



(24) 



;i - RfRe){1 - RfRe 

since phase shift due to GW has to be converted into photo current at the photo detector 
by multiplying a constant factor 24 1. 

Let us summarize this section. We gave detailed description of SRI and identified the 
specific forms of GW signal and noise. The GW signal h^^^^\f) completely corresponds 
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to h{f) in Sec. II, which can be calculated using T in Eq. (I20l) . The power spectrum 
of shot noise can be calculated with Eq. (!23|) . given experimental parameters. In the next 
section, we will investigate the dependence of overlap reduction function 7 on the geometrical 
configuration of two SRIs. 

IV. DEPENDENCE OF SENSITIVITY ON THE RELATIVE LOCATIONS BE- 
TWEEN TWO DETECTORS 

The SNR is significantly influenced by the relative location of two detectors through 
the overlap reduction function when the wavelength of GW is comparable to the size of 
a detector. In this section, we will analyze the overlap reduction function in detail and 
investigate the optimal locations of two detectors in an experiment for the detection of 
GWB at 100 MHz. 

The overlap reduction function can be calculated numerically from Eq. ( fT3|l using the 
arm response function T given in Eq. (1201) . where the phase factor e27n/n Ax/c pig^yg 
important role. To see this, we consider the four configurations of detectors and calculate 
7(/). The results are shown in Fig. [3] with the case of "exact" and "long wavelength limit". 
The former is calculated with the full arm response function T. The latter is calculated 
with T = 1, which is just plotted for reference, though the approximation is not valid 
around 100 MHz. Each configuration of detectors is characterized by the relative position 
AX = Xi — X2 and the relative angle /3. Note that Xj is the position vector of Mi of 
i-th detector (see Fig. [1]). In the case of (a) ideal, 7(/) rapidly decreases even though the 
detectors are completely colocated and coaligned, because we defined the 7(/) including arm 
response functions. Case (b) T-shaped has a behavior similar to that of (a) for the same 
reason. The arm response function is needed to take the effect of the phase change of GW 
at high frequencies into account. Cases (a) and (b) are similar, but have a subtle difference 
since the arms of detectors are at different locations and experience different phases of GW. 
As a result, the overlap of case (b) is a little worse than ). In the cases of (c) crossed 

and (d) stacked, 7(/) also decrease more rapidly than in the long wavelength limit. This 
is because the contribution of the GW phase change at high frequencies is added to that 
in the long wavelength limit. Therefore, we cannot obtain 7(/) = 1 at 100 MHz with the 
detectors where detector size and GW wavelength are comparable. It follows that SNR is 
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worsened by a factor of (0.377)"^ ~ 2.65 at 100 MHz in contrast to the case where long 
wave approximation is vahd, even if the two detectors' configuration is optimal. One may 
expect unit response to be obtained by constructing much smaller detectors. However, the 
total response of detectors 5$ is worsened since the resonant frequency of GW signal also 
depends on the detector size and is shifted upward. Thus, this loss of sensitivity due to the 
phase change of GW is inevitable. 

Next, we will fix the frequency at 100 MHz and consider 7(/). In fact, SRI has a narrow 
frequency band and what we are most interested in is 7 at 100 MHz. In Fig. HJ the location 
of one detector is fixed, while the other detector is located at the same site (AX = 0) and 
the directions of arms are rotated. In this case, the magnitude of 7(/) oscillates, however, 
it has the maximum peak at /5 = and the minimum peak at /5 = vr. It is intuitive that the 
GW signal is better correlated when the directions of arms are coaligned. In Figs. [5] and [6l 
the angle of detectors is fixed and the locations are translated. In an initially coaligned case 
(/3 = 0) in Fig. [5l as expected, 7(/) has its maximum at AX = and keeps the moderate 
value in the range of AX = ±0.2 m. In an initially reversed case [jS = tt) in Fig. [6l an 
interesting feature can be seen. When the detector is translated to the direction (u + v)/v^, 
the peak of 7(/) is shifted. This is because the overlap of the two detectors is better when 
their arms are overlapped geometrically like (c) in Fig. [31 

V. SENSITIVITY TO GWB 

We will describe what the best location of detectors with respect to the sensitivity is, 
and calculate the sensitivity achievable with correlation analysis. From the results obtained 
above, the best location is obviously colocated and coaligned, and gives 7(/)|ioomhz ~ 0.377. 
As shown in Fig. this value is hardly changed in the range of AX = ±0.2 m for coaligned 
detectors. In an experiment, it is impossible to put the detectors in completely colocated 
and coaligned location because of the restricted experimental space of the optics. However, 
experimental detector configuration does not significantly affect 7(/) if the detectors are 
nearly colocated and coaligned. Therefore, we fix it to 7opt = 0.377. 

As for the power spectral density of noise in an SRI, one can calculate from Eq. (!23l) . 
We select the arm length of a detector as L = 0.75 m so that the GW signal resonates at 
100 MHz. With experimental parameters, cu = 1.77 x lO^^rads"^ and 77 = 1, the power 
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FIG. 3: (color online) Overlap reduction function in the case of four detector configurations. 
Each setup is (a) ideal, AX = 0, /? = 0, (b) T-shaped, AX = 0, /3 = 7r/2, (c) crossed, AX = 
(L/2,L/2,0), 13 = TT, (d) stacked, AX = (0,0, L/2), (3 = 0. The "exact" means the calculation 
with arm response function T and the "long wavelength limit" T = 1. The latter is not valid 
around 100 MHz, but merely plotted for comparison. Note that the sign of 7(/) in (b) is inversed 
for convenience of comparison. 



spectral density of noise 3^ around 100 MHz is 

1.60 X 10^^ ^ 



Pi{f) ^ 4.65 X 10 



-42 



IW 



Hz" 



1,2. 



(25) 



"'(/) J \Io 

The factor a' is called the optical amplification factor in a cavity and gives a' ~ 1.6 x 10^ with 
the reflectivity of the recycling mirror, = 0.99996, and the reflectivity of the other three 
mirrors, = (0.99998)^. The bandwidth is ~ 2 kHz with these reflectivities. Substituting 
Pi{f) and 7opt into Eq. (|T71) . and assuming that observation time is T = lyr and that 
^gw(/) has a fiat spectrum around 100 MHz (which is sufficient for practical purposes |29| ) . 
one can calculate the sensitivity of two SRIs to GWB and obtain /iinn^g 



iioo^gw ^ 1.4 X 10 



14 



13 



0.5 




pi/2 pi 3pi/2 2pi 



/3[rad] 

FIG. 4: Overlap reduction function when the detector is initially colocated and coaligned and is 
rotated at the same location. (3 is the rotation angle. 



0.4 




AX [m; 



FIG. 5: (color online) Overlap reduction function when the detector is initially colocated and 
coaligned (/5 = 0) and is translated in certain directions. Each curve means the direction of 
translation. +y, 0) is the direction of (u + v)/ -v/2, (— x, +y, 0) is the direction of (u — v)/\/2, 
and (0,0, z) is the direction perpendicular to the uv plane. 

VI. CONCLUSIONS 

In this paper, we investigated the optimal location of two SRIs whose arms are 0.75 m, and 
derived its sensitivity to GWB in correlation analysis. At 100 MHz, the wavelength of GW 
is comparable to the size of a detector. This means that the GW response of the detector 
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FIG. 6: (color online) Overlap reduction function when the detector is initially colocated and 
reversed (/? = vr) and is translated in certain directions. Each curve means translation in the same 
direction as shown in Fig. [5l 



is less effective than one in the long wavelength limit, and that the location of detectors 
significantly affects the sensitivity. We included the effect due to the finite size of a detector 
into the overlap reduction function and evaluated it. As a result, SNR is worsened by a factor 
of (0.377)^^ ~ 2.65 at 100 MHz in contrast to the case where long-wave approximation is 
valid, even if the two detectors' configuration is optimal. This is because the phase of GW 
changes during light's roundtrip in an arm. SNR also depends on the relative distance and 
angle between detectors. However, SNR is almost optimal value if two detectors are in 
the range of ±0.2 m and coaligned. Using this configuration and experimentally achieved 
parameters of two SRIs, one can achieve the sensitivity to GWB, h^i^^^^ ^ 1.4 x 10^^ 
corresponding to the amplitude h ~ 10~^^Hz~^/^. This constraint on GWB would be much 
tighter than that obtained by current direct observation, h ~ 10~^^Hz~^/^ 1^23,]. However, 
to reach indirect constraint on GWB by big-bang nucleosynthesis, further improvement of 
the sensitivity is required. 
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[31] Note that our derivation of shot noise is slightly different from our previous paper 

we have incorporated the angular response to GW into shot noise. However, here an angular 

averaging is contained in the GW signals. 
[32] Note that this is not the squared strain-amplitude noise ordinarily used in the noise curve of 

ground-based detectors. To convert, one has to multiply Eq. (j25p by (27r/r/ sin27r/r)^, in this 

paper, which is incorporated into the GW signal. 
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